Direct Evidence for Neutrino Flavor Transformation from Neutral- Current 
Interactions in the Sudbury Neutrino Observatory 



Q.R. Ahmad, 17 R.C. Allen, 4 T.C. Andersen, 6 J.D. Anglin, 10 J.C. Barton, 11 - Q E.W. Beier, 12 M. Bercovitch, 10 
J. Bigu, 7 S.D. Biller, 11 R.A. Black, 11 I. Blevis, 5 R.J. Boardman, 11 J. Boger, 3 E. Bonvin, 14 M.G. Boulay, 9 - 14 
M.G. Bowler, 11 T.J. Bowles, 9 S.J. Brice, 9 - 11 M.C. Browne, 17 ' 9 T.V. Bullard, 17 G. Buhler, 4 J. Cameron, 11 
Y.D. Chan, 8 H.H. Chen, 4 >0M. Chen, 14 X. Chen, 8 - 11 B.T. Cleveland, 11 E.T.H. Clifford, 14 J.H.M. Cowan, 7 
D.F. Cowen, 12 G.A. Cox, 17 X. Dai, 11 F. Dalnoki-Veress, 5 W.F. Davidson, 10 P.J. Doe, 17 ' 9 ' 4 G. Doucas, 11 
M.R. Dragowsky, 9 * 8 C.A. Duba, 17 F.A. Duncan, 14 M. Dunford, 12 J.A. Dunmore, 11 E.D. Earle, 14 - 1 S.R. Elliott, 17 - 9 

H.C. Evans, 14 G.T. Ewan, 14 J. Farine, 7 ' 5 H. Fergani, 11 A.P. Ferraris, 11 R.J. Ford, 14 J.A. Formaggio, 17 
M.M. Fowler, 9 K. Frame, 11 E.D. Frank, 12 W. Frati, 12 N. Gagnon, 11 ' 9 - 8 - 17 J.V. Germani, 17 S. Gil, 2 K. Graham, 14 
D.R. Grant, 5 R.L. Hahn, 3 A.L. Hallin, 14 E.D. Hallman, 7 A.S. Hamer, 9 - 14 A. A. Hamian, 17 W.B. Handler, 14 
R.U. Haq, 7 C.K. Hargrove, 5 P.J. Harvey, 14 R. Hazama, 17 K.M. Heeger, 17 W.J. Heintzelman, 12 J. Heise, 2,9 
R.L. Helmer, 16 - 2 J.D. Hepburn, 14 H. Heron, 11 J. Hewett, 7 A. Hime, 9 M. Howe, 17 J.G. Hykawy, 7 M.C.P. Isaac, 8 
P. Jagam, 6 N.A. Jelley, 11 C. Jillings, 14 G. Jonkmans, 7 ' 1 K. Kazkaz, 17 P.T. Keener, 12 J.R. Klein, 12 
A.B. Knox, 11 R.J. Komar, 2 R. Kouzes, 13 T. Kutter, 2 C.C.M. Kyba, 12 J. Law, 6 I.T. Lawson, 6 M. Lay, 11 
H.W. Lee, 14 K.T. Lesko, 8 J.R. Leslie, 14 I. Levine, 5 W. Locke, 11 S. Luoma, 7 J. Lyon, 11 S. Majerus, 11 H.B. Mak, 14 
J. Maneira, 14 J. Manor, 17 A.D. Marino, 8 N. McCauley, 12 ' 11 A.B. McDonald, 14 ' 13 D.S. McDonald, 12 K. McFarlane, 5 
G. McGregor, 11 R. Meijer Drees, 17 C. Mifflin, 5 G.G. Miller, 9 G. Milton, 1 B.A. Moffat, 14 M. Moorhead, 11 
C.W. Nally, 2 M.S. Neubauer, 12 F.M. Newcomer, 12 H.S. Ng, 2 A.J. Noble, 16 ' 5 E.B. Norman, 8 V.M. Novikov, 5 
M. O'Neill, 5 C.E. Okada, 8 R.W. Ollcrhead, 6 M. Omori, 11 J.L. Orrell, 17 S.M. Oser, 12 A.W.P. Poon, 8 ' 17 ' 2 ' 9 
T.J. Radcliffe, 14 A. Roberge, 7 B.C. Robertson, 14 R.G.H. Robertson, 17 ' 9 S.S.E. Rosendahl, 8 J.K. Rowley, 3 
V.L. Rusu, 12 E. Saettler, 7 K.K. Schaffer, 17 M.H. Schwendener, 7 A. Schiilke, 8 H. Seifert, 7 ' 17 ' 9 M. Shatkay, 5 
J.J. Simpson, 6 C.J. Sims, 11 D. Sinclair, 5 ' 16 P. Skensved, 14 A.R. Smith, 8 M.W.E. Smith, 17 T. Spreitzer, 12 
N. Starinsky, 5 T.D. Steiger, 17 R.C. Stokstad, 8 L.C. Stonchill, 17 R.S. Storey, 10 ' | B. Sur, 1 ' 14 R. Tafirout, 7 
N. Tagg, 6 ' 11 N.W. Tanner, 11 R.K. Taplin, 11 M. Thorman, 11 P.M. Thornewell, 11 P.T. Trent, 11 Y.I. Tserkovnyak, 2 
R. Van Berg, 12 R.G. Van de Water, 9 ' 12 C.J. Virtue, 7 C.E. Waltham, 2 J.-X. Wang, 6 D.L. Wark, 15 ' 11 ' 9 
N. West, 11 J.B. Wilhelmy, 9 J.F. Wilkerson, 17 ' 9 J.R. Wilson, 11 P. Wittich, 12 J.M. Wouters, 9 and M. Yeh 3 

(SNO Collaboration) 

1 Atomic Energy of Canada, Limited, Chalk River Laboratories, Chalk River, ON KOJ 1J0, Canada 
2 Department of Physics and Astronomy, University of British Columbia, Vancouver, BC V6T 1Z1 Canada 
3 Chemistry Department, Brookhaven National Laboratory, Upton, NY 11973-5000 
"* Department of Physics, University of California, Irvine, CA 92717 
5 Carleton University, Ottawa, Ontario K1S 5B6 Canada 
6 Physics Department, University of Guelph, Guelph, Ontario NIG 2W1 Canada 
7 Department of Physics and Astronomy, Laurentian University, Sudbury, Ontario P3E 2C6 Canada 
8 Institute for Nuclear and Particle Astrophysics and Nuclear Science 
Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 
9 Los Alamos National Laboratory, Los Alamos, NM 8754-5 
10 National Research Council of Canada, Ottawa, ON K1A 0R6, Canada 
11 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, 0X1 3RH, UK 
Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104-6396 
13 Department of Physics, Princeton University, Princeton, NJ 08544 
lif Department of Physics, Queen's University, Kingston, Ontario K7L 3N6 Canada 

5 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, and 
University of Sussex, Physics and Astronomy Department, Brighton BN1 9QH, UK 
16 TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada 

17 Center for Experimental Nuclear Physics and Astrophysics, and 
Department of Physics, University of Washington, Seattle, WA 98195 
(Dated: 3 February 2008) 

Observations of neutral-current v interactions on deuterium in the Sudbury Neutrino Obser- 
vatory are reported. Using the neutral current, elastic scattering, and charged current reac- 
tions and assuming the standard 8 B shape, the v e component of the 8 B solar flux is <f> e = 
1.76l°;o5(stat.)+° °g (syst.) x 10 6 cm~ 2 s" 1 for a kinetic energy threshold of 5 MeV. The non-v e 
component is (fr^ — 3.41+q' 4 ! (stat.)^' 48 (syst.) x 10 6 cm~ 2 s -1 , 5.3a greater than zero, providing 
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strong evidence for solar v a flavor transformation. The total flux measured with the NC reaction is 
0nc = 5.091q 43 (stat.)^Q *g (syst.) x 10 6 cm _2 s _1 , consistent with solar models. 

PACS numbers: 26.65.+t, 14.60.Pq, 13.15.+g, 95.85.Ry 



The Sudbury Neutrino Observatory (SNO) detects 8 B 
solar neutrinos through the reactions: 

v e + d -> p+p + e" (CC), 
v x + d -> p + n + u x (NC), 
v x +c~^ v x + e" (ES). 

The charged current reaction (CC) is sensitive exclusively 
to electron- type neutrinos, while the neutral current re- 
action (NC) is equally sensitive to all active neutrino fla- 
vors [x — e, /i,r). The elastic scattering reaction (ES) 
is sensitive to all flavors as well, but with reduced sensi- 
tivity to and v T . Sensitivity to these three reactions 
allows SNO to determine the electron and non-electron 
active neutrino components of the solar flux jl| . The CC 
and ES reaction results have recently been presented . 
This Letter presents the first NC results and updated CC 
and ES results from SNO. 

SNO H is a water Cherenkov detector located at a 
depth of 6010 m of water equivalent in the INCO, Ltd. 
Creighton mine near Sudbury, Ontario, Canada. The de- 
tector uses ultra-pure heavy water contained in a trans- 
parent acrylic spherical shell 12 m in diameter to de- 
tect solar neutrinos. Cherenkov photons generated in the 
heavy water are detected by 9456 photomultiplier tubes 
(PMTs) mounted on a stainless steel geodesic sphere 17.8 
m in diameter. The geodesic sphere is immersed in ultra- 
pure light water to provide shielding from radioactivity 
in both the PMT array and the cavity rock. 

The data reported here were recorded between Nov. 2, 
1999 and May 28, 2001 and represent a total of 306.4 live 
days, spanning the entire first phase of the experiment, in 
which only D2O was present in the sensitive volume. The 
analysis procedure was similar to that described in 
PMT times and hit patterns were used to reconstruct 
event vertices and directions and to assign to each event 
a most probable kinetic energy, T c g. The total flux of 
active 8 B solar neutrinos with energies greater than 2.2 
MeV (the NC reaction threshold) was measured with the 
NC signal (Cherenkov photons resulting from the 6.25 
MeV 7 ray from neutron capture on deuterium.) The 
analysis threshold was T c g> 5 MeV, providing sensitivity 
to neutrons from the NC reaction. Above this energy 
threshold, there were contributions from CC events in 
the D2O, ES events in the D2O and H2O, capture of 
neutrons (both from the NC reaction and backgrounds), 
and low energy Cherenkov background events. 

A fiducial volume was defined to only accept events 
which had reconstructed vertices within 550 cm from the 
detector center to reduce external backgrounds and sys- 
tematic uncertainties associated with optics and event 



reconstruction near the acrylic vessel. The neutron re- 
sponse and systematic uncertainty was calibrated with 
a 252 Cf source. The deduced efficiency for neutron cap- 
tures on deuterium is 29.9 ± 1.1% for a uniform source 
of neutrons in the D2O. The neutron detection effi- 
ciency within the fiducial volume and above the en- 
ergy threshold is 14.4%. The energy calibration was 
updated from Q with the 16 N calibration source 
data and Monte Carlo calculations. The energy response 
for electrons, updated for the lower analysis threshold, 
was characterized as a Gaussian function with resolution 
a T = -0.0684 + 0.331^ + 0.0425T e , where T e is the 
true electron kinetic energy in MeV. The energy scale 
uncertainty is 1.2%. 

The primary backgrounds to the NC signal are due to 
low levels of uranium and thorium decay chain daughters 
( 214 Bi and 208 T1) in the detector materials. These activ- 
ities generate free neutrons in the D2O, from deuteron 
photodisintegration (pd), and low energy Cherenkov 
events. Ex-situ assays and in-situ analysis of the low 
energy (4 — 4.5 MeV) Cherenkov signal region provide 
independent uranium and thorium photodisintegration 
background measurements. 

Two ex situ assay techniques were employed to deter- 
mine average levels of uranium and thorium in water. Ra- 
dium ions were directly extracted from the water onto ei- 
ther MnO x or hydrous Ti oxide (HTiO) ion exchange me- 
dia. Radon daughters in the U and Th chains were subse- 
quently released, identified by a spectroscopy, or the ra- 
dium was concentrated and the number of decay daughter 
(3-a coincidences determined. Typical assays circulated 
approximately 400 tonnes of water through the extrac- 
tion media. These techniques provide isotopic identifica- 
tion of the decay daughters and contamination levels in 
the assayed water volumes, presented in Fig. |l] (a). Sec- 
ular equilibrium in the U decay chain was broken by the 
ingress of long-lived (3.8 day half-life) 222 Rn in the exper- 
iment. Measurements of this background were made by 
periodically extracting and cryogenically concentrating 
222 Rn from water degassers. Radon from several tonne 
assays was subsequently counted in ZnS(Ag) scintillation 
cells @j. The Radon results are presented (as mass frac- 
tions in g(U)/g(D 2 0)) in Fig. 0(b). 

Independent measurements of U and Th decay chains 
were made by analyzing Cherenkov light produced by 
the radioactive decays. The f3 and (3-j decays from the 
U and Th chains dominate the low energy monitoring 
window. Events in this window monitor 7 rays that pro- 
duce photodisintegration in these chains (E 7 > 2.2 MeV). 
Cherenkov events fitted within 450 cm from the detector 
center and extracted from the neutrino data set provide 
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a time-integrated measure of these backgrounds over the 
same time period and within the fiducial volume of the 
neutrino analysis. Statistical separation of in situ Tl and 
Bi events was obtained by analyzing the Cherenkov sig- 
nal isotropy. Tl decays always result in a (5 and a 2.614 
MeV 7, while in this energy window Bi decays are domi- 
nated by decays with only a f3, and produce, on average, 
more anisotropic hit patterns. 

Results from the ex situ and in situ methods are con- 
sistent with each other as shown on the right hand side 
of Figs. 0(a) and 0(b). For the 232 Th chain, the weighted 
mean (including additional sampling systematic uncer- 
tainty) of the two determinations was used for the anal- 
ysis. The 238 U chain activity is dominated by Rn ingress 
which is highly time dependent. Therefore the in-situ de- 
termination was used for this activity as it provides the 
appropriate time weighting. The average rate of back- 
ground neutron production from activities in the D2O 
region is 1.0 ± 0.2 neutrons per day, leading to 44tg de- 
tected background events. The production rate from ex- 
ternal activities is OT^g neutrons per day, which leads 
to 27 ±8 background events since the neutron capture ef- 
ficiency is reduced for neutrons born near the heavy wa- 
ter boundary. The total photodisintegration background 
corresponds to approximately 12% of the number of NC 
neutrons predicted by the standard solar model from 8 B 
neutrinos. 

Low energy backgrounds from Cherenkov events in the 
signal region were evaluated by using acrylic encapsu- 
lated sources of U and Th deployed throughout the de- 
tector volume and by Monte Carlo calculations. Prob- 
ability density functions (pdfs) in reconstructed vertex 
radius derived from U and Th calibration data were 
used to determine the number of background Cherenkov 
events from external regions which either entered or mis- 
reconstructed into the fiducial volume. Cherenkov event 
backgrounds from activities in the D2O were evaluated 
with Monte Carlo calculations. 

Table | shows the number of photodisintegration and 
Cherenkov background events (including systematic un- 
certainties) due to activity in the D2O (internal region), 
acrylic vessel (AV), H2O (external region), and PMT ar- 
ray. Other sources of free neutrons in the D2O region 
are cosmic ray events and atmospheric neutrinos. To 
reduce these backgrounds, an additional neutron back- 
ground cut imposed a 250-ms deadtime (in software) fol- 
lowing every event in which the total number of PMTs 
which registered a hit was greater than 60. The num- 
ber of remaining NC atmospheric neutrino events and 
background events generated by sub-Cherenkov thresh- 
old muons is estimated to be small, as shown in Table |. 

The data recorded during the pure D 2 detector phase 
are shown in Figure [| These data have been analyzed 
using the same data reduction described in 0, with the 
addition of the new neutron background cut, yielding 
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FIG. 1: Thorium (a) and uranium (b) backgrounds (equiv- 
alent equilibrium concentrations) in the D2O deduced by in 
situ and ex situ techniques. The MnO x and HTiO radio- 
chemical assay results, the Rn assay results, and the in situ 
Cherenkov signal determination of the backgrounds are pre- 
sented for the period of this analysis on the left-hand side of 
frames (a) and (b). The right-hand side shows time-integrated 
averages including an additional sampling systematic uncer- 
tainty for the ex situ measurement. 



TABLE I: Neutron and Cherenkov background events. 



Source 


Events 


D2O photodisintegration 




H2O + AV photodisintegration 


27± 8 8 


Atmospheric i/'s and 




sub-Cherenkov threshold /x's 


4± 1 


Fission 


< 1 


2 H(q, a)pn 


2 ±0.4 


17 0(a,n) 


< 1 


Terrestrial and reactor P's 


1±? 


External neutrons 


< 1 


Total neutron background 


78 ± 12 


D2O Cherenkov 


20 + H 


H2O Cherenkov 


0+4 

J -3 


AV Cherenkov 


6t 3 6 


PMT Cherenkov 


16 + ^ 


Total Cherenkov background 


451- 



4 



TABLE II: Systematic uncertainties on fluxes. The experi- 
mental uncertainty for ES (not shown) is -4.8,-1-5.0 percent, 
f denotes CC vs NC anti-correlation. 



Energy scale f -4. 2, +4. 3 -6. 2, +6.1 -10. 4, +10.3 

Energy resolution f -0.9, +0.0 -0.0, +4. 4 -0.0, +6. 8 

Energy non-linearity t ±0.1 ±0.4 ±0.6 

Vertex resolution f ±0.0 ±0.1 ±0.2 

Vertex accuracy -2. 8, +2. 9 ±1.8 ±1.4 

Angular resolution -0.2, +0.2 -0.3, +0.3 -0.3, +0.3 

Internal source pd f ±0.0 -1.5, + 1.6 -2.0, +2. 2 

External source pd ±0.1 -1.0, + 1.0 ±1.4 

D 2 Cherenkov f -0.1, +0.2 -2. 6, +1.2 -3. 7, + 1.7 

H 2 Cherenkov ±0.0 -0.2, +0.4 -0.2, +0.6 

AV Cherenkov ±0.0 -0.2, +0.2 -0.3, +0.3 

PMT Cherenkov f ±0.1 -2.1, +1.6 -3.0, +2. 2 

Neutron capture ±0.0 -4.0, +3. 6 -5. 8, +5. 2 

Cut acceptance -0.2, +0.4 -0.2, +0.4 -0.2, +0.4 



Experimental 


X 


icertainty 


-5. 2, +5. 2 


-8. 5, +9.1 


-13. 2, + 14.1 


Cross section 






±1.8 


±1.3 


±1.4 



2928 events in the energy region selected for analysis, 5 
to 20 MeV. Fig. |||(a) shows the distribution of selected 
events in the cosine of the angle between the Cherenkov 
event direction and the direction from the sun (cos 9q ) 
for the analysis threshold of T e ff> 5 MeV and fiducial 
volume selection of R < 550 cm, where R is the recon- 
structed event radius. Fig. |||(b) shows the distribution of 
events in the volume-weighted radial variable (i?/i?Av) 3 j 
where i?AV = 600 cm is the radius of the acrylic ves- 
sel. Figure ||(c) shows the kinetic energy spectrum of the 
selected events. 

In order to test the null hypothesis, the assumption 
that there are only electron neutrinos in the solar neu- 
trino flux, the data are resolved into contributions from 
CC, ES, and NC events above threshold using pdfs in T e g, 
cos#0, and (R/Rav) 3 , derived from Monte Carlo calcu- 
lations generated assuming no flavor transformation and 
the standard 8 B spectral shape JsJ . Background event 
pdfs are included in the analysis with fixed amplitudes 
determined by the background calibration. The extended 
maximum likelihood method used in the signal decompo- 
sition yields 1967.7±|o % CC events, 263.6t;|; 4 ES events, 
and 576.5^4gg NC events Jl2"| , where only statistical un- 
certainties are given. Systematic uncertainties on fluxes 
derived by repeating the signal decomposition with per- 
turbed pdfs (constrained by calibration data) are shown 
in Table |Tj. 

Normalized to the integrated rates above the kinetic 
energy threshold of T c ff> 5 MeV, the flux of 8 B neutri- 
nos measured with each reaction in SNO, assuming the 
standard spectrum shape [[j) is (all fluxes are presented 



in units of 10 6 cm 2 s 1 ): 
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Electron neutrino cross sections are used to calculate all 
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FIG. 2: (a) Distribution of cos 6 Q for R < 550 cm. (b) Dis- 
tribution of the volume weighted radial variable (R/ Rav) 3 ■ 
(c) Kinetic energy for R < 550 cm. Also shown are the 
Monte Carlo predictions for CC, ES and NC + bkgd neutron 
events scaled to the fit results, and the calculated spectrum 
of Cherenkov background (Bkgd) events. The dashed lines 
represent the summed components, and the bands show ± la- 
uncertainties. All distributions are for events with T e g>5 
MeV. 
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fluxes. The CC and ES results reported here are consis- 
tent with the earlier SNO results @ for T cff >6.75 MeV. 
The excess of the NC flux over the CC and ES fluxes 
implies neutrino flavor transformations. 

A simple change of variables resolves the data di- 
rectly into electron (<j> e ) and non-electron {4>^ T ) compo- 
nents Q , 

</> e = 1.76±ro5(stat.)±°;°09 (syst.) 
V = 3.41±^(stat.)+^ (syst.) 

assuming the standard 8 B shape. Combining the sta- 
tistical and systematic uncertainties in quadrature, 0^ 
is 3.4llg g4, which is 5.3er above zero, providing strong 
evidence for flavor transformation consistent with neu- 
trino oscillations ^ [|. Adding the Super-Kamiokande 
ES measurement of the 8 B flux (l0| 0|f = 2.32 ± 
0.03(stat.)+°;°7 (syst.) as an additional constraint, we 
find (f)^iT = 3.45lQ g2, which is 5.5a above zero. Fig- 
ure ^ shows the flux of non-electron flavor active neutri- 
nos vs the flux of electron neutrinos deduced from the 
SNO data. The three bands represent the one standard 
deviation measurements of the CC, ES, and NC rates. 
The error ellipses represent the 68%, 95%, and 99% joint 
probability contours for tp e and </> Mr . 

Removing the constraint that the solar neutrino energy 
spectrum is undistorted, the signal decomposition is re- 
peated using only the cosOq and (R/Rav) 3 information. 
The total flux of active 8 B neutrinos measured with the 
NC reaction is 

<f$P = 6.42tl;g(stat.)^:l (syst.) 

which is in agreement with the shape constrained value 
above and with the standard solar model prediction |Tl| ] 
for 8 B, SSM = 5.05^. 

In summary, the results presented here are the first di- 
rect measurement of the total flux of active 8 B neutrinos 
arriving from the sun and provide strong evidence for 
neutrino flavor transformation. The CC and ES reaction 
rates are consistent with the earlier results j| and with 
the NC reaction rate under the hypothesis of flavor trans- 
formation. The total flux of 8 B neutrinos measured with 
the NC reaction is in agreement with the SSM prediction. 

This research was supported by: Canada: NSERC, In- 
dustry Canada, NRC, Northern Ontario Heritage Fund 
Corporation, Inco, AECL, Ontario Power Generation; 
US: Dept. of Energy; UK: PPARC. We thank the SNO 
technical staff for their strong contributions. 
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